Pair annihilation in laser pulses: optical vs. XFEL regimes 
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We discuss the theory and phenomenology of pair annihilation, within an ultra-short laser pulse, to a single 
photon. The signature of this process is the uni-directional emission of single photons with a fixed energy. We 
show that the cross section is significantly larger than for two-photon pair annihilation in vacuum, with XFEL 
parameters admitting a much clearer signal than optical beams. 



Introduction. Prospects for probing the quantum vacuum 
using laser light at facilities such as ELI and the European 
XFEL cuiTently attract a great deal of attention. Taking into 
account both the finite duration and spatial extent of modern 
laser pulses is challenging, but good analytic progress can be 
made if one assumes the laser is not too tightly focussed; it 
may then be modelled by a time dependent electric field, or 
a plane wave. Scattering in such backgrounds is very well 
understood, especially when the pulse is treated as having 
an effectively infinite duration [h-iH- These results have re- 
cently been re-evaluated, following the advent of ultra-short 
laser pulses, in order to expose the effects of finite pulse ge- 
ometry on laser-particle scattering 13^81. One may therefore 
still ask if, within this well known topic, there are unique ex- 
perimental signatures from which one could easily identify 
strong-field effects. 

We discuss such a signal here, based on electron-positron 
annihilation. Unlike in vacuum, pair annihilation in a laser 
pulse may produce a single photon, both the energy and mo- 
mentum of which are fixed. Photons from each annihilation 
event are emitted in the same direction and with the same en- 
ergy. In other words, repeated (idealised) events would them- 
selves produce a laser-like beam of photons. Since the pro- 
duced photons carry half the energy of the incoming, accel- 
erated, particles, they have energies well above the electron 
rest mass: our annihilation process would produce a gamma 
ray beam! Even if gamma-ray lasers are still a pipe dream, 
we note that the most promising route to their developement 
is via the annihilation of molecular positronium, as has been 
demonstrated in the lab |9|. The process we investigate here 
is essentially an 'accelerator based' version of that idea. 

We begin by summarising the properties of our process, be- 
fore giving its exact cross section in a laser pulse of finite du- 
ration, using strong-field QED. This is then investigated for a 
variety of realistic optical and XFEL beam parameters. For 
related recent investigations see flOl FTl I . 

Kinematics. In vacuum, momentum conservation implies 
that an electron-positron pair must annihilate to (at least) two 
photons lfT2l . In a laser field, though, the one-photon chan- 
nel opens. This process is rather novel since there is only one 
allowed final state, modulo polarisations: momentum conser- 
vation alone determines the scattering products, as we now 
describe. We model the laser by a plane wave function of 
k.x, with the laser momentum lightlike. A pair with mo- 
menta and p'^ (mass m) annihilate, using up some quantity 
Ik^ of momentum from the laser fields, and produce a sin- 




FIG. 1: Left: annihilation to two photons in vacuum. Right: the 
one-photon channel in a laser pulse. Double lines indicate dressing 
of the fermions by the laser field. 

gle photon of momentum A:^ . Momentum conservation gives 
Pn + p'^ + Ik^ ~ k'^ . There are four components to this equa- 
tion, but a real photon has only three (momentum) degrees of 
freedom. This implies that the energy drawn from the laser is 
fixed: indeed, one finds after squaring up that I is determined 
by the incoming momenta as / = — (m^ + p.p')/k.{p + p'). 
This constraint on I gives the right counting of degrees of free- 
dom. That I is negative may be interpreted as energy being 
given up to the laser Note that I is not an integer, in gen- 
eral. (It effectively becomes so only for periodic backgrounds, 
after explicitly subtracting that part of I which generates the 
electron mass shift 11131114 :1.) To illustrate, let the pulse move 
down the z axis, so kfi = mv( 1,0,0, 1) with v = (o/m the di- 
mensionless laser frequency. The pair are introduced from the 
transverse directions with equal energy E = my and opposite 
momenta. One easily finds from the above that ^ = — 7/v and 
k!^ = 1,0, 0,-1). Hence, half of the incoming energy is 
given up to the field and half (always) produces a backscat- 
tered photon with high frequency o' = my. 

The cross section. The laser field strength is F^vik.x) — 
fi{k.x)(ki^a[, — a'^kv). With A^ as given above, the transverse 
polarisation vectors are a'^ = (moo / e) di^ for / = 1,2. Our nor- 
malisation is such that the invariant, dimensionless amplitude 
flo is the peak laser intensity. The two functions fi (k.x) and 
/2 (k.x) give the shape of the laser pulse. (The derivative on 
fj is for convenience, since it is the fj themselves which ap- 
pear in the cross section.) Define /a = ff +/! ■ We take the 
fields to vanish outside of k.x G [Q,2nN], so that 2nN is the 
Lorentz invariant pulse duration and the number of laser 
cycles. Introducing the pulse envelope function g{k.x), such 
that g vanishes smoothly at the pulse edges and \g\ < 1, we 
parameterise fi {x) = g{x) cos(x) and f2{x) = g{x) sin(x). 

In vacuum, a scattering cross section is given by the proba- 
bility for a process to occur per unit volume per unit time (i.e. 
one divides out the interaction volume, which is the volume 
of spacetime, after squaring the 5-matrix element), divided by 
the incoming particle flux. The cross section a for our process 



is defined analogously: the only difference is that the interac- 
tion volume is the pulse volume. To proceed, we introduce 
functions of the external momenta: 
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along with u = k.k'^ /{4-k.pk.p') and = InP'a^u/k.k' , which 
appear in three 'Bessel-hke' integrals (/ — 1,2,3), 




FIG. 2: Collision geometry in the lab. The pair travel anti-parallel to 
the laser (blue, momentum fc^) with an offset angle B. The emitted 
photon (red, k'^) is backscattered along the laser direction. 
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Repeated indices are summed over j E {1,2,3}. Define also 
Bo = {(^jBj)/! 1 15 1 and the combination 

Z = 2|Bo|^ +«o(2« - l)(2|Bi \^+2\B2\^-BoBl-BlBi) . 

The flux is / = V {p-p'Y ^ ™d cross section is 
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Some brief remarks on known cases. In a periodic laser field 
(A? = oo), ([3]) is a delta comb (T1I21: to make sense of this one 
must add various effects of finite duration by hand lfT3l . This 
problem does not appear at all if one includes the physical 
pulse duration from the outset, as we have done. In a peri- 
odic field, ([3]) is a function of the effective electron mass, but 
this plays no explicit role here, as has been observed for other 
processes in short pulses QIlll: the energy of the produced 
photon, for example, does not depend on the effective mass. 

In a constant (crossed) field, (|3]l vanishes. To understand 
why, compare with another laser process in the constant field 
limit, say nonlinear Compton scattering Q. In that process, 
the interaction volume diverges, which appears to suppress the 
rate. However, it was shown in [IJ that the integrals over the 
final states also contribute an infinite and equivalent factor, so 
the total cross section remains finite. On the other hand, the 
final state space in our process is trivial and independent of the 
details of the background, so no such cancellation can occur in 
general. Hence the emission rate vanishes. We now examine 
a for various beam and particle parameters. 

Low intensity optical beams. We begin by following |fT6l , 
which considered an optical laser with intensity oq = 1 • • - 5. 
This is outside the perturbative regime (oq < 1) but low by 
modern standards: ELI beams will have aq ^ 10"*. With- 
out resorting to perturbation theory, |16 | provided an elegant 
method for evaluating the functions (|2]i. The integrals therein 
seem suited to evaluation by stationary phase, but the kine- 
matics of the process forbid the exponents from having any 
stationary phase points. One may instead separate the inte- 
grand into slowly and rapidly varying parts if the number of 
laser cycles is large, and then sets the rapid oscillation scale 
ifTTl . (The method is particularly suited to the XFEL regime, 
as we will shortly see.) This leads to an expansion of the in- 
tegrand of Q in Bessel functions J„, similarly to the case of 



a periodic plane wave, for details see ifTSl §101]. The remain- 
ing rapidly oscillating exponent now admits stationary phase 
points, which are the solutions of 
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For given and p'^ (i.e. given I and aa) there is a finite 
range of negative n which contributes. For each n the station- 
ary phase points are calculated, and a sum over n performed 
to obtain the Z?,. (We have verified this approach using the 
powerful numerical methods of HPil.) The rapid decay of the 
Bessel functions /„ with increasing n means that the most sig- 
nificant contributions to a come from low n. Since > 
implies l — n>0, such contributions approximately require 
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where j3 — -y/ 1 — 1 / as usual. The leading factor in ^ in- 
creases with decreasing laser frequency, and so very small an- 
gles 9 are required to keep the value of / low. An appreciable 
signal will therefore require the setup shown in Fig. |2] with 
the pair's momenta a small angle 9 from being anti-parallel to 
the laser direction ll20l . Assuming an optical frequency (0 — 5 
eV (v = 10^^), intensity ^ 1 • • • 5 and a large number of cy- 
cles N ^ I (this does not model a modern pulse, see below), 
along with high energy pairs, y ^ 10^, (7 exceeds that for the 
two-photon channel in vacuum by over an order of magni- 
tude, for a range of collision angles 9 |fT3l,IT6|. This is where 
we encounter the difficulty in observation of our process: the 
required angle is 0° < 10^^ because of the low optical fre- 
quency. Such precision seems experimentally challenging, at 
best, so we turn to alternative laser parameters. We are nat- 
urally led to consider a high frequency XFEL beam. We will 
see in this regime that not only does the cross section continue 
to dominate over that in vacuum, but that one gains several or- 
ders of magnitude in "experimental ease". 

The XFEL regime. We assume a 100 fs pulse at X-ray fre- 
quency (x> ~ ^ 10^ eV (v = 1 /lOO), generated by the European 
XFEL at DESY |2ll|22, and 7 = 80 (already available at the 
DESY FZD hnac). This pulse contains A? = 1.2 x 10^ cycles 
due to its high frequency: the 'slowly varying phase' approxi- 
mation introduced above is therefore extremely well suited to 
the XFEL regime, as coiTections go like 1 /N. To be concrete 
we choose the envelope function g{x) = sm'^{x/2N) following 
|5|. Our initial intensity is oq = 1/10, which will shortly be 
increased. 
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FIG. 3: Cross section a x s {h = c = I) in a 100 fs XFEL pulse of 
intensity ao = I/IO, v — 1/100, with y = 80. Top: the cross section 
dominates over that in vacuum (red, bottom line) for a narrow angu- 
lar range (vertical scale curtailed for easier comparison). Bottom: a 
closer view of the above peak. CT exhibits a rich substructure. 



Two views of the resulting cross section, multiplied by 
Mandelstam s= {p + p'Y, are shown in Fig.jsjas a function of 
the collision angle Q. Since our intensity is low we compare 
the cross section with that for two-photon annihilation in vac- 
uum. There are clear quantitative and qualitative differences: 
while the two photon process is supported over all 0, the one- 
photon channel is very strongly peaked over a narrow angular 
range. In the lower panel of Fig.|3]we 'zoom in' on the cross 
section to expose its oscillatory substructure, which is an ef- 
fect of finite pulse duration as discussed in |4, 7, 14|. While 
the support of the peak is now offset from = by roughly 
half a degree, its angular range remains very narrow. 

The visibility of the signal is significantly improved by rais- 
ing the XFEL intensity to the edge of the perturbative domain, 
flo = 1 (other parameters as above). The result is shown in the 
top panel of Fig. |4] the increased intensity leads to an appre- 
ciable signal over an angular range of ^ 1°. The appearance 
of disjoint signals (each with a rapidly oscillating substruc- 
ture) corresponds to (|4]) admitting solutions for — n — 1 ... 3 
over the given Q range: increasing the intensity allows more 
n to contribute since the denominator in (|4ji goes like Aq. The 
precise form of a is sensitive to the intensity, as illustrated in 
the lower panel of Fig. |4]where we consider the 'goal' XFEL 
intensity ao = 10, and observe a signal over ^ 2°. Here, many 
n give overlapping contributions to a 1 13|. Similar intensity- 
dependent behaviour is seen in nonlinear Compton scattering 
||23I . The two-photon result is shown in Fig.|4]for comparison, 
but other processes may contribute in the nonperturbative do- 



main flo > 1, for example the two-photon channel in the laser 
(which is suppressed by a factor of flg perturbatively). The lit- 
erature does not seem to contain a complete treatment of this 
process, though the case of a periodic field may be extracted 
by crossing symmetry from |24|, and prospects for coherent 
X-ray generation via this process are discussed in |25|. 

Despite the challenge in detecting the signals presented 
here, we emphasise that the cross section is significantly 
greater than that in vacuum for the shown 9 ranges. Further- 
more, going to XFEL beams improves the visibility of the sig- 
nal by four orders of magnitude: the angular range in Fig.|4]is 
measured in whole degrees, 10^ times larger than that in the 
optical regime. 

High intensity optical beams. The next generation of opti- 
cal lasers will have extremely high intensities of flo ^ 10^, and 
ultrashort duration ~ ^(1)- As such, neither a perturbative 
expansion in the intensity flo nor the above "large A^" method 
are applicable. Numerically, we have not found an apprecia- 
ble signal for the geometry shown in Fig. |2] and reasonable 
collision angle. We therefore return to the transverse collision 
discussed in the introduction, which is a more natural experi- 
mental setup. We take flo = 10"* and O) = 1 eV (v = 2 x 10"^), 
modelling an ELI strength optical laser With this, each of the 
terms in the exponents of Q is large: 



-10" 7, 



as 



lo'Vr- (6) 



Assuming high energy incoming pairs, say y ^ 10"*, each of 
the parameters in ^ is around ^ =M. The functions ^ 




FIG. 4: At higher intensity the cross section is supported over a 
wider angular range of ~ 2°. Top: at ao = 1 we again see a rapidly 
oscillating substructure. Bottom: at uq = 10 many n must be summed 
over to obtain a. Omitting the oscillating substructure, we have 
drawn the envelope functions of the first six, dominant, contributions. 
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in this regime therefore take the form 
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Since the full exponent contains no stationary phase points 
(h' ^ 0), one may generate an asymptotic expansion of these 
integrals using the standard 'integration by parts' method: 
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This is just the statement that the functions (j2]i are dominated, 
for short pulses and at high energy, high intensity, by edge 
effects. If the pulse turns on and off symmetrically then = 
at the pulse edges, and the above term is zero. Going to higher 
orders in the expansion, one encounters the derivatives of /, at 
the pulse edges. At least the first derivative must vanish if we 
are to have a smooth pulse. The first nonvanishing term in the 
expansion ([s]) therefore goes like M^^^ ' , with r the order of 
the first non-vanishing derivative of /, at the pulse edge: the 
smoother the pule, the smaller the resulting cross section. 

Relevance to cascade formation. It has been suggested 
that a cascade of particles could be triggered even by a single 
pair creation event in an intense laser pulse Il26lj28l . Cascade 
codes use cross sections calculated in a constant, crossed field, 
as this is how any laser field looks locally, at sufficiently high 
intensity ||T]. As discussed above, the cross section Q in such 
a field vanishes fT '^l, and so it seems that one-photon annihi- 
lation is not particularly relevant to cascades. Even with finite 
pulse duration, we have seen that CT is appreciable only for 
particular collision geometries with small angular tolerances, 
and such fine tuning is of course impossible in a cascade event. 
(Density effects ||29l and higher order processes such as the 
two-photon channel f2A\ are a different matter, and must be 
discussed elsewhere.) 

Conclusions. We have discussed pair annihilation to one 
photon, within a laser pulse. This process has a trivial final 
state space and cannot occur in vacuum. We have provided a 
modern treatment which incorporates, from the outset, physi- 
cal effects due to finite pulse duration. This removes the theo- 
retical ambiguities which appear for periodic fields. We have 
analysed the cross section for beam parameters which are, or 
will soon be, realisable at the European XFEL and ELI. 

We have seen that the cross section for our process can be 
significantly larger than that of the two-photon channel in vac- 
uum, but that its observation requires fine tuning of the colli- 
sion geometry. The best situation we have found is for a rea- 
sonably intense XFEL beam, in which the angular precision 
required is measured in degrees. This is an improvement from 
thousandths of degrees in the optical regime and is due to the 
high XFEL frequency. After submitting this paper the preprint 
l30 1 appeared which supports these conclusions using a quan- 
tum kinetic approach. 
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